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We propose an implementation of the quantum fast Fourier transform algorithm in an entangled 
system of multilevel àtoms. The Fourier transform occurs naturally in the unitary time evolution of 
energy eigenstates and is used to define an alternate wave-packet basis for quantum information in the 
atom. A change of basis from energy levels to wave packets amounts to a discrete quantum Fourier 
transform within each atom. The algorithm then reduces to a series of conditional phase transforms 
between two entangled àtoms in mixed energy and wave-packet bases. We show how to implement such 
transforms using wave-packet control of the internal states of the ions in the linear ion-trap scheme for 
quantum computing. 
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The discrete quantum Fourier transform, 
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which links two sets of states each labeled by integers, 
occurs in many applications of quantum computing [Q. 
It is central to Shor's algorithm for prime factorization 
[2|, which has applications in públic- key cryptography 
Coppersmith [Q describes an efficient algorithm for 
implcmenting this transform whcn N is a power of 2, 
achieving an exponential speed-up over the classical fast 
Fourier transform (FFT) algorithm ||. 

Advances in quantum FFT technology have been mo- 
tivated by the difficulty of implcmenting quantum lògic 
between macroscopically distinct two-level systems, or 
qubits. The difficulty arises from decoherence, the loss 
of coherence in a quantum superposition due to coupling 
with the environment. It is known that one-qubit gates 
alone, interspersed by classical measurements, suffice to 
build the quantum FFT |(| , but this result is hard to im- 
plement in practice. Approximate simulations of DFTjy 
have also been proposed Q,|t| , and are known to be more 
tolerant to phase fluctuations in the two-qubit gates when 
applied in the context of Shor's algorithm ||. 

In this paper, we consider an analogue of the exact 
quantum FFT algorithm based on multi-valued quan- 
tum lògic ||, and propose a novel realization of DFTjv 
in multilevel atòmic systems using wave-packet con- 
trol methods. The advantage of using d > 2 comput- 
ational levels in each atom is that the number of 
àtoms needed for the algorithm is reduced by a factor 
of log 2 d. For example, d — 8 levels stores three qubits 
of information in each atom, requiring only Q/3 àtoms 
for computing DFT^ for N = 2 Q . Since fewer àtoms 
are needed, the multilevel approach minimizes the 



decoherence associated with the macroscopic entangle- 
ment of these àtoms, and enables a seale-up in the im- 
plementation of the quantum FFT. 

In section ^ we show that the elementary operations 
needed for the algorithm are a Fourier transform of the 
d levels in each atom, DFT^, and a phase gate that 
couples two àtoms together. The d-level Fourier trans- 
form takes the place of the Walsh-Hadamaard transform, 
which plays a prominent role in binary quantum compu- 
tation pcj| . The phase gate involves a conditional cou- 
pling between two entangled àtoms, and is more suscepti- 
ble to decoherence in implementation. As the number of 
phase gates in the quantum FFT seales as the square of 
the number of àtoms, the reduction in the latter in a mul- 
tilevel implementation is advantageous from a coherence- 
time standpoint. 

We propose to implement DFT^ in the atom by a 
change of computational basis, as described in section |^. 
The Fourier transform occurs naturally in quantum me- 
chanics in relating complementary representations and 
we show that this can be useful for computational pur- 
poses. A dual Fourier basis for atòmic energy levels con- 
sists of localised electron wave packets at discrete Urnes 
in one Kepler orbit about the nucleus |yj. Individual 
elements in the wave-packet basis can be addressed by 
short làser pulses that interact with the electron when it 
is near the atòmic core. A change of basis from energy 
levels to wave packets effectively accomplishes DFT^ in 
the atom. 

The quantum FFT then reduces to a sequence of con- 
trolled phase gates between two àtoms in hybrid bases, 
evolving the phases of wave-packet states in one atom 
conditional on energy levels in the other. In section 
we consider a method for implementing such a gate in the 
linear ion-trap quantum lògic scheme proposed by Cirac 
and Zoller |12[. A multilevel phase-gate protocol in this 
scheme involves a sequence of làser pulses applied to two 
ions in the trap. 
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2. MULTI- VALUED QUANTUM FFT 

In a system of Q = log 2 N qubits, DFTjv can be cons- 
tructed using only two kinds of binary gates [§],[Ï3| . These 
are the single-qubit Walsh-Hadamaard transform 
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acting on qubit to, and the two-qubit controlled phase 
gate 
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acting on qubits l and to, where (j) — Ti/2 m ~ l . Specifically, 
it can be shown that except for a reversal of bits in the 
final output, 

DFTjv = (Aq_iBq_ 2 ,Q-i)(^Q-2-Bq_3 i q_i_Bq_3 ! Q_2) 

. . . (AiB ,q_i5o,q-2 ■ ■ • B ,i)A , (4) 

where the sequence of gates on the right side is applicd 
from left to right. The total number of gates is Q(Q + 
l)/2 = 0[Q 2 }, so this is an efficient process. 

To describe the multi-valued quantum FFT, we gener- 
alize the gates A m and Bi m to multilevel systems. Each 
d-level system is referred to as a qudit. The states \a) 
and \c) in Eq. (g) can be written as a tensor product of 
q = \og d N qudits, 



\a) = \a q -x,a q -2, a ), 
a m = 0, 1, . . . , d — 1 for all m, 
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and similar ly for |c). The numbers a m represent the dig- 
its of a in base d. The number of qudits q in the tensor 
product is less than the number of qubits Q by a factor 
of log 2 d, 
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which reduces the number of àtoms needed for the 
algorithm. The multi-valued analogue of the Walsh- 
Hadamaard transform A m is a d-level Fourier transform, 
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which mixes the d states in the mth qudit, |0),|1), 
. . . , \ d— 1), with phases determined by the Fourier kernel. 
The phase gate Bi m generalizes to the two-qudit gate 
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which is a diagonal transformation that advances the 
phase of each of the d 2 states in the two-qudit basis, 
|0, 0), |0, 1), . . . , \d— 1, d— 1), by an amount determined 
by the vàlues of both qudits, a; and b m . For d = 2 and 
ai, b m = or 1 in Eqs. (M) and (^]), we recover the binary 
gates in Eqs. (@) and (|j) respectively. 

Given the above definitions for A m and £>; m , we show 
that a sequence of gates similar to that in Eq. (^) simu- 
lates DFTjv on an g-qudit register. In the multi-valued 
case, for N = d q , 

DFTyv = (^,_lB g -2, g -l)(^,-2B,-3,,-lB,-3,,-2) 

...(4iflo, a -iflo, g -2.·.flo,iMo· (9) 

Based on an argument given by Shor Q for the binary 
quantum FFT, we consider the matrix element of DFTyv 
between two arbitrary states \a) and |c), 
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and show that the sequence of gates in Eq. (Q) has the 
same matrix element as above, but between states |o) 
and |6), where \b) is defined as the 'dit-reversed' version 
of \c), 



\b) = \b q -!,b q -2, ■ ■ -,b ), 
= \c ,Ci, . . . ,Cç_l>. 
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The least signihcant place in b becomes the most signifi- 
cant place in c, and vice versa. A reversal of qudits can 
be performed efficiently using multi-valued permutation 
gates H , or else we can simply read out the final state in 
the reverse order. 

To detcrmine the amplitude Aé 1 ^ of going from 
|a,_i,a ? _2,...,ao) to \bq-i, b q s, . . . , b ) under the se- 
quence of gates in Eq. (||), consider each set of gates 
separated parenthetically in this sequence. First A m 
transforms \a m ) to \b m ) in the mth qudit with ampli- 
tude (1/yd) exp(i2Tta m b m /d). This is followed by all 
the gates £>; m for to > l, each of which adds a phase 
2iraib m /d m ~ l+1 to \ai,b m ) without mixing states. The 
net modulus A of the transition amplitude between |o) 
to \b) is thus determined by the product of the A m gates, 
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The net phase cj> can be separated into two parts, that 
due to the A m gates and that due to the Bi m gates, 
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Since the first term amounts to setting l = to in the sec- 
ond term, we can combine the terms by replacing to > l 
with to > l in the second summation. From Eq. (fill) , 
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FIG. 1. Illustration of quantum FFT gates (based on Ref. 
[13]). The shaded squares represent the qudits transformed 
by Am, implemented as a change from the energy-level to the 
wave-packet basis. 

we have b m = c g _i_ m . Defining m' = q — 1 — to, the 
summation over m > l becomes one over m! < q — l, 
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Including m! > q — l tcrms in the second summation 
above will not afïect the phase since these extra terms 
are integer múltiples of 2-7T. Hence, the two summations 
decouple to give 
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where we have used N = d q and identified a and c in 
their base-d notation. From Eqs. ( |Ï2| ) and (|ï^), we see 
that the net amplitude of going from \a) to \b) under 
the sequence of gates in Eq. (JsJ) is identical to that of 
going from \a) to |c) under DFTat. Thus, to within a 
reversal of qudits between |ò) and |c), the q(q+ï)/2 gates 
in Eq. (^|) simulates a quantum Fourier transform on a 
g-qudit register. 

A graphical illustration of the quantum FFT is shown 
in figurc [!]. The first three passes through the algorithm 
(to = q — 1, q— 2, q— 3) are shown, corresponding to the 
first three sets of gates in Eq. (|9|). During each pass, an 
A m gate first mixes the d states in qudit to, illustrated 
by a shading of the respective square in the figure, fol- 
lowed by a sequence of Bi m gates that couple all shaded 
squares with the next unshaded one. In each pass, the 
A m gate enables a ci-point Fourier transform that is re- 
peated efHciently d m times by the conditional Bi m gates, 
achieving exponential speed-up over the classical FFT. 
The speed-up is made possible by the tensor product na- 
ture of quantum entanglement |Ï4j . In the multi- valued 
case, this corresponds to a d-ary tree decomposition of 
DFTjy in terms of unitary operations. 

At the complction of the quantum FFT, each qudit is 
read out by measurement. For a single multilevel sys- 
tcm, this yields one value per qudit, corresponding to 
one of the d computational levels. Information stored in 
a superposition of these levels is lost upon measurement, 



analogously to the situation in a two-level system. The 
measured output of the multi-valued algorithm is thus 
always a product state of the q qudits, corresponding to 
a classical number in base d. 



3. ATÒMIC FOURIER TRANSFORM 

Consider the implementation of the transform A m in 
an atom with d computational levels. As shown in 
Eq. (0), this transform uniformly mixes all the levels in 
the atom with phases determined by the Fourier kernel. 
In the basis of energy levels, this requires precise con- 
trol of the relative phases of the levels. Such control 
is not feasible for large numbers of levels in the energy 
basis. However, we can regard this as a problem in wave- 
packet control. Atòmic wave packets are superpositions 
of energy levels with different phase relations among the 
levels. We proposo to implcmcnt DFT^ in the atom by 
means of a dual computational basis composed of wave 
packets pjj. 



3.1. Energy and wave packet bases 

The Fourier transform occurs naturally in the change 
of state representation from coordinate to momentum in 
quantum mechanics. Although time is not an observable, 
we can also speak of an uncertainty relation between en- 
ergy and time. In this case, the Fourier kernel appears 
in the unitary time-evolution operator, which relates the 
continuous dynamics of the bound atòmic state to its dis- 
crete energy spectrum. A larger number of energy levels 
in the superposition leads to a greater localization in the 
wave packet. By appropriately discretizing the dynam- 
ics, we can define a wave-packet basis in the atom that 
is related to the energy-level basis by a discrete quantum 
Fourier transform. 

Consider radial wave packets Pq| , which are superpo- 
sitions of Rydberg energy levels with low angular mo- 
mentum. These levels have long radiative life times, ap- 
proaching a millisecond for principal quantum number 
n > 100. Take d energy levels ccntcrcd at n with angular 
momentum l = 1 to represent the computational basis in 
the atom, 
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\j) v = \n + j,l,0), 
-d/2 + 1, -d/2 + 2,..., d/2, 



(16) 



where the subscript v denotes a state in the energy-level 
basis. We have assumed above that n is an integer and 
d is an even number for simplicity, but the arguments 
below are easily extended to non-integer n and odd d. 
A uniform superposition of the d levels corresponds to a 
radially-localized wave packet in space whose time evo- 
lution is given by 
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where hujj is the energy of the jth level in the superposi- 
tion, and Hujo is the mean energy of the wave packet cor- 
responding to principal quantum number n. To separate 
the classical and revival dynamics of the wave packet, we 
expand u)j — luq in a Taylor series in j = n — fi, 



Lüj — LÜQ = 2lT 
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The Kepler period Tk = 2irn 3 (in atòmic units) mea- 
sures the round-trip time for the wave packet traveling 
between the inner and outer turning points of the classi- 
cal or bit. This corresponds to a radial shell of probability 
distribution varying periodically in size. The revival time 
T rev and the super-revival time T SI describe higher-order 
effects such as dispersion and revivals. 

Define an orthogonal basis of wave-packet states corre- 
sponding to d discrete times during the classical Kepler 
evolution of the radial wave packet Q, 



l fc )r = -^ÇexpH27rjÀ:/d)|j) v 
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where k = —d/2 + 1, —d/2 + 2, . . . , d/2, and the subscript 
r denotes a state in the wave-packet basis. Thus k = 
corresponds to the wave-packet state centered at the in- 
ner turning point of the classical orbit, and k = í\k\ 
correspond to wave-packet states moving in opposite di- 
rections at some intermediate location in the orbit, as 
illustrated in figure ^. 

The energy-level basis \j) v is related to the wave-packet 
basis |fc) by a <i-level Fourier transform, 



A m \j) u = ^=^exp(z2^7d)|j') ly 
y/d -, 

= \k= -j) T , 
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where we have used Eqs. @j and ©. This suggests 
that A m can be implemented in the atom by a change 
of computational basis from that of energy eigenstates to 
that of wave-packet states. A change of basis does not 
involve any free time evolution or real-time processing in 
the atom, which means that there is no computational 
cost to realizing the A m gates in this manner. 

To understand the effect of the basis change on the 
algorithm, we refer again to figure 0. Recali that the 
shading of each square in the figure corresponds to 
the application of a A m gate to that qudit in the ar- 
ray, and the following Bi m gates in that pass through 
the algorithm couple all the shaded squares with the 
next unshaded onc. If A m is regar ded as a change 
of basis from energy levels to wave packets, then the 
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FIG. 2. The Fourier conjugate of the energy-level basis \ j) 
consists of wave-packet states \k) evenly distributed in time 
around a classical orbit. A radial wave-packet basis can be 
thought of as an ensemble of such orbits with different orien- 
tations for the ellipses. 

shaded squares are to be read in the wave-packet basis. 
Consequently, each Bi m gate involves different bases for 
the two qudits l and m in the transformation. The quan- 
tum FFT thus reduces to a series of conditional two-qudit 
phase transforms Bi m performed in hybrid bases. 

In the atòmic case, Bi m can be regarded as phase 
shifts applied to each wave-packet state |fc) in the mth 
atom conditional on each energy eigenstate \j) in the 
Ith atom. In view of Eq. (|20|) , we rewrite Eq. m) as 
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We describe a protocol for implementing Bi m in a linear 
ion trap in section ^. This requires coherent control of the 
wave-packet basis in the target atom, which we discuss 
below. 



3.2. Coherent wave packet control 

Consider an arbitrary Rydberg state in the wave- 
packet basis, 



|^(t)) =exp(-twoí)5^6*(*)|fc) 7 



(22) 



where bk are slowly varying amplitudes from which we 
have removed the average frequency ujq corresponding to 
the mean Rydberg level fi in Eq. ([Ï6|). Since the wave- 
packet states are not stationary, the amplitudes evolve 
in time. During a Kepler period, the periòdic motion of 
the radial wave packet corresponds to a cyclic permuta- 
tion in the amplitudes, 
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where m is an integer and k — m is taken modulo d. At 
later times, higher-order terms in the Taylor expansion 
of Eq. ( |ïlj| ) become significant, and the dispersion of the 
Rydberg wave function mixes the amplitudes b k nontriv- 
ially. However at the revival times T Tev , the wave func- 
tion reforms into the original state and nearly recovers 
the initial distribution of amplitudes in the wave-packet 
basis. 

An applied làser field interacts strongly with a Rydberg 
wave packet when it is near the atòmic core. This phe- 
nomenon underlies the excitation and photo-ionization 
of radial wave packets |ïr| , and we use this as a means 
for coherent control of individual amplitudes b k in the 
wave-packet basis. The idea is to use short làser pulses 
to transfer a chosen amplitude to the ground state for 
time-resolved processing. 

Consider a broadband làser pulse with a spectral width 
~ d/Tx that couples all the Rydberg levels in Eq. JÏ6| ) 
to the ground state \g) in the atom. If the pulse is 
transform-limited, it has a temporal width less than 
Tk / d, and we can to good approximation ignore the free 
Kepler evolution of the wave-packet amplitudes shown 
in Eq. (p3|). Only the amplitude òo changes significantly 
during the pulse, corresponding to the wave-packet state 
|0) r nearest the atòmic core. This state undergoes Rabi 
oscillations with the ground state pd[ , 



f(t) f2 exp(-ïA í)ò , 



6 o — 2 /(í)^oexp(iA í)ò 9 . 
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where f(t) is the pulse profile, Ao = ujq — lo is the de- 
tuning of the center frequeney u> of the pulse from the 
average Rydberg frequeney, and ílo is proportional to 
the average of the Rabi freqüències íl g j for the ground- 
to-Rydberg transitions, 



(26) 



A localized wave packet behaves classically for t ~ Tr, 
and the strong coupling to the làser field near the core 
can be understood as a large momentum transfer to the 
electron near the nucleus. The ratc at which energy is 
absorbed from the field E is proportional to p-E, and the 
electron momentum p is maximum at the inner turning 
point. 

The two-level system of Eqs. ( f2~i|) and (^5|) allows selec- 
tive wave-packet processing in the atom. In particular, 
it allows phase control of individual wave-packet states 
in the Rydberg basis, as needed to implement B[ m . To 
see this, note that a 7r-pulse of duration less than Txjd 
de-excites only that part of the Rydberg wave function 
associated with a single wave-packet amplitude b k , creat- 
ing a 'dark' wave packet in its place in the Rydberg basis. 
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FIG. 3. Phase-gate implementation. (a) Ion trap with 
pulses applied to ions l (control) and m (target). (b) Level 
scheme in each ion, with Rydberg levels j and auxiliary levels 
g and e. Three types of pulses used: Ui(u)j — v x ,'ïï), V m (^o, i"), 
and U m (üü e — v x — A, 2tt). See text for details. 



The phase of bk in the ground state can be controlled 
by a narrow-band 27r-pulse that couples this state to an 
auxiliary energy level in the atom (level e in figure |^). 
The phasc-shiftcd amplitude can be restored to the Ry- 
dberg basis at a time commensurate with when the dark 
wave-packet state returns to the atòmic core. This ac- 
complishes a phase shift of a selected wave-packet state 
in the target atom. 

Since the wave-packet basis is not stationary, the time 
between successive interactions with the Rydberg basis is 
dictated by the free atòmic time seales. This time interval 
has to an integer number of Kepler periods during which 
the Rydberg state has not incurred much dispersion, or 
alternately, when the free time evolution undergoes a re- 
vival in the Schròdinger picture. This ensures that the 
decomposition of the Rydberg wave function in the wave- 
packet basis has not changed appreciably between pulses. 



4. PHASE GATE IN THE LINEAR ION TRAP 

Consider an implementation of the two-qudit phase 
gate Bi m in the linear ion-trap scheme for quantum 
computing |p^| . Assuming that d energy levels in each 
trapped ion represent a qudit, consider a series of làser 
pulses applied to ions / and m in the trap, as illustrated 
in figure 0(a). Our goal is to control the phases of the 
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wave-packet states \k) T in the mth ion conditional on 
the energy eigenstates \j) in the Ith ion, as required by 
Eq. ©• 

The ions are assumed to be in the vibrational ground 
state and oscillate synchronously in the center-of-mass 
normal mode in the trap. Assuming that the interaction 
ficld has a standing-wave pattern along the trap axis, two 
kinds of interactions have been proposed in this scheme, 
labeled U and V ]Ï2fl . The V interaction arises when 
the ion is at the antinode of the standing wave, and the 
làser resonantly couples two internal states in the ion 
according to the unitary evolution operator 



V(í) = exp 



(27) 



where Í2 is the Rabi frequency and à is the lowering op- 
erator for the atòmic transition. Alternately, the U inter- 
action arises when the ion is at the node of the standing 
wave and the làser is detuned off resonance to an atòmic 
transition by the trap frequency v x . We consider the low- 
est two trap states, |0) and For atòmic levels g and e, 
where e is the upper level, we find that the states |<jf)|l) 
and |e)|0) are coupled by the unitary operator 



U(t) = exp 



-it [<j'a 

V? 2 
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where à is the trap lowering operator and q is the num- 
ber of ions in the trap. The Lamb-Dicke parameter rj is 
defined as 



V = k x 



2mv x ' 
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where m is the mass of each ion, and k x is the wave vector 
along the trap axis. The unitary evolution in Eq. ( pq ) is 
vàlid in the limit that r\ <C 1. 

Consider the two-ion Rydberg wave function at some 
time íq when the trap has been initialized to |0), 
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where we use the abbreviation \j') v \k') = \j',k'), and 
the summations over j' and k' run over the d components 
of the energy-level and wave-packet bases in ions l and 
m respectively. The coefhcients Cjty are the Schròdinger 
picture amplitudes whose free time evolution has two 
contributions, one due to the phase evolution of the en- 
ergy levels in the Zth ion, and another due to the periòdic 
evolution of the wave-packet amplitudes in the mth ion. 
We have to kcep thesc contributions in mind as we pursue 
a phase gate in the hybrid basis. 

When the fcth wave-packet element in the mth ion 
is near the atòmic core, the methods described in sec- 
tion [3 3.2 can be used to transfer the corresponding am- 
plitudes c.j'k to the ground state \g) in the ion. This is 
done by applying a broadband 7r-pulse of the V type, de- 
noted by V m (üJo,Tr) in figure ||(b). The pulse spectrum 



is centered on the mean frequency ujq and has a dura- 
tion less than Tx/d that is an integer múltiple of n/ílo- 
This only affects the wave-packet state nearest the atòmic 
core, denoted by k' — k, and leaves the two ions in the 
state 
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The second term in Eq. (|3lj) represents that part of the 
wave function in the mth ion that is left in the Rydberg 
manifold, and we suppress this term briefly. The first 
term corresponds to the fcth wave-packet state that has 
been de-excited, which can be written as 



|* fc (íi)) = E c ^(*i)IÍ'=5)|0) 



= Cj*(ti)|i,5>|0) + E cfk(h)\f,g)\0), (32) 

where a particular energy level j is taken out of the f 
summation. We seek to de-excite this level to the ground 
state in the Zth ion, conditional on exciting the trap. This 
is done by applying a narrow-band 7r-pulse of the U type 
to the Ith ion, which has a pulse duration that is an 
integer múltiple of 7r / '(r)Q g j / \/q). The làser frequency 
is tuncd to ojj for the jth Rydberg level. This pulse 
is labeled Ui(uij —v x ,ii) in figure |3|(b), and transforms 
|*fc(*i)> to 

\* k (t 2 ))=c jk (t 2 )\g ig )\l) + J2cj'k(t2)\f,g)\0), (33) 

where the coefficients have evolved in phase from t\ to t 2 
due to the free time evolution of the energy levels in the 
Schròdinger picture. Equation ( |33| ) shows that the trap 
is excited only when both ions are in the ground state 
\g, g), corresponding to the initial state \j,k) at time ío- 
Hence, the U pulse has created entanglement between 
the trap state and the internal states of the two ions. 

To implemcnt Bi m , we need to shift the phase of | j, k) 
by 4>jk according to Eq. (|2~ï|). In state |\&fc(Í2)), this cor- 
responds to evolving the phase of \g, g)\l) by (f>jk without 
affecting the other basis states. To do this, consider the 
auxiliary level e in the mth ion shown in figure 0(b) . Ap- 
plying a U pulse of 2n duration couples states \g)\l) and 
|e)|0) in the mth ion. For a làser detuning of A, the in- 
teraction phase of \g, g)\l) evolves by an integer múltiple 
of 7r(l + A/Qg e ), which can be controlled to achieve <j>jk- 
This pulse is denoted U m (u> e — i/ x — A, 2n) in the figure, 
and transforms lí'fc^)) to 



|*fc(*a)) = c jk (t 3 ) exp(z^ fc ) \g, g)\l) 

+ J2 c Mh)\f,g)\o), 



(34) 



G 



giving a controlled phase shift 4>jk to the state \g, g) |1) as 
desired. We now reverse the operation that took us from 
|*k(íi)) to |*fc(í 2 )) by applying Ui (u)j - v x ,tt) again to 
the Zth ion, creating 

|*fc(Í4)) = c jk (n)exp(i(j) jk )\j,g)\0} 

+ J2 c fk(t4)\f,9)\0). (35) 

Lastly, the mth ion state \g) is restored to \k) by ap- 
plying V m (lüo , 7r) again at a time that is commensurate 
with when the 'dark' radial wave-packet element corre- 
sponding to fc returns to the atòmic core. Since this is a 
V pulse, it does not affect the trap. The resulting state 
is 

|**(*s)> = c ifc (í 5 )exp(# ife )|j,fc)|0) 

+ J2°Mt5)\j',k)\0). (36) 

The k' 7^ k terms in Eq. ( |3~ï|) are unaffected by the com- 
bination of the five pulses used above. Including their 
contribution to the final wave function, we get 

|*(t 6 ))=c ifc (Í5)exp(i^ fc )|i,fc)|0) 

+ EE c ^(t 5 )\j',k')\Q). (37) 

Comparing Eqs. (|3C]) and (|37|), we see that the sequence 
of five pulses, 

V m (uo,ir) Ui(wj-i/ x ,Tr) U m (u; e - i/ x - A, 2w) 

Ui(üJj-v x ,Tr)V m (üJ ,ir), (38) 

accomplishes a controlled phase shift of the hybrid state 
\j,k) = \j) ll \k) T , as desired. This procedure is repeated 
for each of the d 2 states in the two ions, leading to the 
phase gate Bi m in mixed energy-level and wave-packet 
bases. 

The coefficients Cjiy (£5) are different from Cfk' (£0) due 
to the phase evolution of the energy levels in both ions 
during the time interval £5 — £0. This is also responsible 
for the non-stationarity of the wave-packet basis, which 
makes this procedure sensitive to when the fcth wave- 
packet amplitude in the mth ion is de-excited from, and 
excited to, the Rydberg manifold. The time interval be- 
tween the two V m pulses in the sequence is determined 
by the free atòmic time seales governing the classical or 
revival dynamics of the wave packets. By appropriate 
timing of these pulses, we can implement a conditional 
phasc-shift between the two àtoms. 

5. CONCLUSION 

This paper shows that the quantum FFT can be sim- 
plificd using a multilevel wave-packet approach to its im- 
plementation. In atòmic systems, the bàsic lògic gate is 



a multilevel Fourier transform A m in each atom, which 
is equivalent to a change of basis from energy levels to 
wave packets. Such an implementation takes advantage 
of the natural Fourier transform relation between energy 
and time in quantum mechanics. The FFT then reduces 
to a series of two-qudit phase gates Bi m in hybrid bases, 
which we have considered in the context of the linear 
ion-trap scheme for quantum computing. 

The advantage of the multilevel approach is a reduc- 
tion in the number of entangled quantum systems (e.g. 
trapped ions) by a factor of log 2 d compared to the bi- 
nary case. For the same reason, the number of lògic 
gates needed to simulate DFTjv is fewer by a factor of 
(log 2 d) 2 , as seen by comparing Eqs. (^) and ([}]). How- 
ever, this comes at the cost of larger elementary gates, 
up to (£ 2 -dimensional in the case of Bi m . The trade-off in 
computation time depends very much on the particular 
implementation scheme used, which dictates the physical 
time taken to perform each gate. 

The bottlcneck for the gate operation time in the linear 
ion-trap scheme is the trap frequeney v x [ [171 , typically 
kHz to MHz, which limits the speed of the narrow-band U 
pulses. This is due to the need for selective entanglement 
between the internal state of the ion and its motional 
state in the trap. By comparison, the Rydberg atòmic 
time seales are much faster, typically in the ns to /is 
range, allowing much faster execution times in principio 
for the phase-gate protocol given in section |J 

There are two advantages to a wave-packet approach 
to multilevel processing. One is the quantum Fourier 
transform itself, which is integral to the approach and key 
to quantum computing. The other is the feasibility of the 
control scheme for atòmic systems. Universal control of 
múltiple energy levels in the atom requires múltiple làsers 
tuned to the neighbouring transitions ||. This is much 
easier in the time domain, where wave-packet transforms 
can be achieved by controlling the timing and durations 
of a sequence of làser pulses. The implementation of the 
quantum FFT in multilevel systems is thus made more 
feasible using wave-packet methods. 

This work was supported by the Army Research Office 
through the MURI Center for Quantum Information. 



Electronic Address: amuthuk@optics.rochester.edu 

[1] R. Jozsa, Proc. Roy. Soc. Lond. A 454, 323 (1998). 

[2] P. W. Shor, SIAM J. Comp. 26, 1484 (1997). 

[3] R. L. Rivest, A. Shamir and L. Adleman, Commun. of 

A.C.M. 21, 120 (1978). 
[4] D. Coppersmith, IBM Res. Report RC 19642, T.J. Wat- 

son Research Center, Yorktown Heights, NY, 1994. 
[5] J. W. Cooley and J. W. Tukey, Math. Computation 19, 

297 (1965). 

[6] R. B. Griffiths and C. Niu, Phys. Rev. Lett. 76, 3228 
(1996). 

[7] R. Cleve and J. Watrous, in Proc. l^lst Annual Symp. 



7 



Found. Com v. Sci. (IEEE, Los Alamitos, CA, 2000), p 



[8] 
[9] 

[10] 

nu 



526. Also at yiant-ph/0006004 



A. Barenco et al., Phys. Rev. A 54, 139 (1996). 

A. Muthukrishnan and C. R. Stroud, Jr., Phys. Rev. A 

62, 52309 (2000). 

Bowden, C. M., Chen, G., Diao, Z., and Klappenec- 



ker, A., quant-ph/0007122 at xxx.lanl.gov. 

Muthukrishnan, A., and Stroud, Jr., C. R., |quant- 



[12] 
[13] 
[14] 



[15] 

[16] 
[17] 



ph/0106165 at xxx.lanl.gov. 
J. I. Cirac and P. Zoller, Phys. Rev. Lett. 74, 4091 (1995). 
A. Ekert and R. Josza, Rev. Mod. Phys. 68, 733 (1996). 
R. Jozsa, to appear in Geomètric Issues in the Founda- 
tions of Scienc e, eds. S. Huggett et al. (Oxford Univ. 
Press). Also at |quant-ph/9707034 . 

J. Parker and C. R. Stroud, Jr., Phys. Rev. Lett. 56, 716 
(1986). 

A. ten Wolde et al., Phys. Rev. Lett. 61, 2099 (1988). 
D. Jonathan, M. B. Plenio and P. L. Knight, Phys. Rev. 
A 62, 42307 (2000). 



8 



